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ABSTRACT 

In 1982, Monique and Frangois Spite discovered that the ^Li abundance in the atmosphere 
of old metal-poor dwarf stars in the galactic halo was independent of metallicity and temper- 
ature. Since then, ^Li abundance in the Universe has become a subject of intrigue, because 
there is less of it in Population 11 dwarf stars (by a factor of 3) than standard big bang nucle- 
osynthesis predicts; a discrepancy which is still far from being solved. The most challenging 
features of the plateau is the lack of variability in the -2.8 < [Fe/H] < -2.0 metallicity range 
and the "meltdown" (a drop and an increase in variability) for [Fe/H] < -2.8. Here we show 
how quark-novae (QNe) occurring in the wake of Pop 111 stars, (a few days to a few weeks 
following the core-collapse SN explosions of 20-40Mo progenitors) can elegantly produce an 
A(Li) ~ 2.2 Lithium plateau in Pop 11 (low-mass) stars formed in the pristine cloud swept up 
by the mixed SN+QN ejecta. We also find an increase in the scatter as well as an eventual 
drop in A{Li) below the Spite plateau values for very low metallicity ([Fe/H] < -3) in excel- 
lent agreement with observations. We propose a solution to the discrepancy between the Big 
Bang Nucleosynthesis ^Li abundance and the Spite plateau and list some implications and 
predictions of our model. 

Key words: cosmology: theory - early universe - nuclear reactions, nucleosynthesis, abun- 
dances - supernovae: general - stars: carbon - stars: Population 11, Population 111 



1 INTRODUCTION 

One of the most important discoveries in the study of the chemical 
composition of stars was made in 1982 by Monique and Francois 
Spite, who found an essentially constant ^Li abundance in warm 
metal-poor stars in the galactic Halo (Spite&Spite 1982). The 
lithium abundance shows a plateau with a near constant valuajof 
A(Li) ~ 2.2 for metallicity between [Fe/H]= -2 and -2.8 (and effec- 
tive stellar temperature 5900 K < Tgff < 6300 K). The star-to-star 
scatter in ^Li abundance in these stars is extremely small, of the 
order of ~ 0.05 dex (e.g. Asplund et al. 2006). 

Given the low nuclear binding energy of ^Li and because its 
cross section to proton capture is quite large, it needs only be ex- 
posed to relatively low stellar temperatures (T ^ 2.5 x 10 K) to 
suifer substantial destruction over the many-billion-year life span 
of a halo star. These temperatures are never reached in the stellar 
atmospheres but if the convection zone is deep, ^Li in the atmo- 
sphere is brought to hot deep layers where it is slowly destroyed. 
In warm {T^fi ^ 5900K) metal-poor stars of Spite&Spite (1982) 
the convection is not so deep, and lithium should be preserved"! 



Therefore, the lithium abundance in the atmosphere of these stars 
is expected to reflect the abundance of ^Li in the primitive Galactic 
matter; these stars hold the memory of the chemical composition 
of the interstellar medium or cloud from which they formed. It was 
proposed to identify the ^Li abundance of the plateau to the abun- 
dance produced by the big bang (primordial lithium); i.e. as a relic 
of primordial nucleosynthesi^J However, several theoretical and 
observational results have cast doubt on use of the observed ^Li 
plateau abundance as the primordial value despite its thinness (e.g. 
Ryan et al. 1999; Cyburt et al. 2008). 

Thanks to high precision WMAP measurements of of the cos- 
mological constants (e.g. Dunkley et al. 2009; see also Larson 
et al. 2011; Komatsu et al. 2011), the Big Bang nucleosynthesis 
(BBN) of the elements (^He, "^He, D, ^Li and ''Li ) during the 
first few minutes in the evolution of the Universe, is now known 
with good precision. All elements occur in amounts predicted, ex- 
cept for ^Li. The corresponding primordial ^Li abundance derived 
is A(Li) = 2.72+[5J55 (Cyburt et al. 2008; Coc & Vangioni 2010), 
which is a factor of 3 times higher than the ^Li abundance inferred 
from warm low-mass halo stars. 



* email:rouyed@ucalgary.ca 

' We adopt the standard definition [A/B] = log (Na INb)- log (Na INb)o and 
A{Li) = logiNij^jINH) + 12 where the subscript O refers to solar values. 
^ In standard stellar models, the expected depletion of 'Li for metal-poor 
tumoff stars is negligible, i.e., less than 0.02 dex (Deliyannis et al. 1990; 
Pinsonneault et al. 1992). 



' The mean value of the lithium abundance measured in the old galactic 
stars is also very close to the value measured in u) Cen cluster of A{Li) = 
2.19±0.14 (see Bonifacio et al., 2010); presumably Ci> Cen is the remnant 
of a captured dwarf galaxy. This suggests that A(Li) = 2.1-2.2 could be a 
"universal value of stellar 'Li abundance at low metallicity. 
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There are several discussions in the hterature about how to 
solve this discrepancy between the measured Spite value and the 
WMAP/BBN results (e.g. Fields 2011 for a recent review). Ideas 
involve: (i) a first generation of primordial, zero-metallicity (Pop 
III), stars that has processed (by nuclear burning) and efficiently 
depleted lithium in a substantial fraction of the early halo bary- 
onic matter before low-mass stars on the Spite plateau formed 
(Piau et al. 2006). This pre-galactic ^Li processing idea faces 
some challenges; it does not appear to be able to explain the en- 
tire WMAP/Spite plateau gap (e.g. Sbordone et al. 2010). Fur- 
thermore, it appeals to Pop III stars with masses in the range 10- 
4OM0 which are known to produce an amount of CNO elements 
much higher than the observed abundances in metal-poor halo stars 
(Prantzos 2007); (ii) Other theoretical studies involve non-standard 
BBN models exploring the frontiers of new physics (e.g. Jedamzik 
et al. 2006; Coc et al. 2009; locco et al. 2009; Jedamzik & Pospelov 
2009; Kohri & Santoso 2009). These modifications to the BBN 
consider the decay of unstable supersymmetric particles (e.g. Cy- 
burt et al. 2010) to lower the primordial prediction. For e.g. if the 
late decaying particle is a gravitino, then the Spite plateau and the 
baryonic density determined by WMAP may be reconciled (e.g. 
Jedamzik et al. 2006). Non-standard BBN cannot explain the de- 
pendence of A(LO on metallicity below [Fe/H] — 2.8; (iii) Studies 
involving modifications to the reaction cross-sections for the ^Li 
production during BBN (e.g. Chakraborty, Fields & Olive 2011), 
although large errors in the cross section of the reaction leading 
to ''Li are unlikely (e.g. Cyburt et al. 2004; Angulo et al. 2005). 
Ideas calling for alterations in ^Li production reaction rates in nu- 
clear physics and in particle physics in BBN have been ruled out 
as a cause (e.g. Coc et al. 2010); (iv) Alternatively, the discrepancy 
can in principle be eliminated by assuming that some degree of ^Li 
depletion has occurred. Diffusion causes ^Li (and other elements) 
to sink below the photosphere where it is destroyed at T ^ 2 x 10^ 
K by reactions with protons (e.g., Michaud et al. 1984; Salaris & 
Weiss 2001). The degree and efficiency of the photospheric deple- 
tion of ^Li can be "controlled" by including mixing in diffusion 
models (e.g. Lambert 2004). For example, mixing due to Gravity 
waves (Talon et al. 2002; Charbonnel & Talon 2005), to Rotation 
(e.g., Vauclair 1988; Chaboyer & Demarque 1994; Pinsonneault 
et al. 2002; Talon & Charbonnel 2004) and, turbulence (Richard 
et al. 2005; Piau 2008). Richard et al. (2005) include a specific 
parametrization of turbulent mixing, and managed to produce a flat 
^Li abundance along the plateau, with A(Li) depleted by about 0.4 
dex with respect to the initial value (see fig. 8 in Richard et al. 
2005; see also Mucciarelli et al. 2011). The turbulence is included 
in these models without postulating any physical mechanism re- 
sponsible for it and, its tuning leads to unwanted effects on the 
abundances of elements other than ^Li (e.g. §6.4 in Bonifacio et 
al. 2007). 



1.1 The challenge 

The lithium problem is far from being solved. The puzzle remains 
the same and in some sense we are still facing severe problems in 
our understanding of how lithium was synthesized in the Big Bang 
and subsequently destroyed in stars. It is still unclear whether the 
disagreement results from lack of understanding of stellar lithium 
destruction mechanisms or from non-standard physics in BBN. 
Any new model must face two key questions (see discussion in 
Spite et al. 2012 and references in Fields 201 1): 



• The very low scatter: Based on data of the highest quality, 
there is very low scatter in the ^Li abundances among stars on the 
^Li plateau; the actual scatter may even be smaller than the ob- 
served one (e.g. Asplund et al. 2006). The negligible scatter (^ 0.05 
dex), and apparent lack of slope in the Spite plateau are observa- 
tional constraints, which ^Li depletion models have trouble repro- 
ducing. 

• The meltdown: The meltdown of the Spite plateau is estab- 
lished, but its cause is unclear (Sbordone et al. 2010). Specifically, 

below [Fe/H] 2.8, more and more stars appear to exhibit ^Li 

abundances below the plateau level, while scattering increasing. 
The challenge is to envision a single depletion phenomenon pro- 
ducing a thin, metallicity independent plateau above [Fe/H] — 2.8, 
and a highly scattered, metallicity dependent distribution below it. 



The ^Li problem has become increasingly acute and this "lin- 
gering" situation leaves room for new ideas to be explored. Here 
we appeal to primordial Quark-Novae (hereafter QNe) going-off 
in the wake of Pop III stars (prior to the era of Pop II stars). The 
core-collapse (Type-II) SN explosion of a Pop III star with mass 
in the 20Mq < Mprog. < 40Mo range would leave behind a neu- 
tron star massive enough to experience an explosive transition to a 
quark star; the QN (Ouyed et al. 2002; Vogt et al. 2004; Ouyed et 
al. 2005; Keranen et al. 2005; Niebergal et al. 2010). A QN explo- 
sion ejects the neutron-rich outer layers of the parent neutron star. 
This ejecta expands away from the quark star (the compact rem- 
nant) at relativistic speeds with an average Lorentz factor Fqn - 10; 
i.e. a typical neutron energy iiQN = 10 GeV (Keranen et al. 2005; 
Ouyed et al. 2005). The interaction between the neutron-rich QN 
ejecta and the preceding (Pop III) type-II SN ejecta leads to spal- 
lation products with abundances reflective of extremely metal-poor 
stars (Ouyed 2013). Spallation (i.e. destruction) of the ^^Ni in the 
ejecta of the (Pop III) type II SN by the QN neutrons leads to iron- 
impoverishment in the cloud swept up by the combined SN+QN 
ejecta. The Pop II stars formed in these clouds will appear metal- 
poor with [Fe/H] decreasing with fdelay (the time delay between 
the SN explosion and the subsequent QN explosion; see Figure 1 
in Ouyed 2013). In this paper, we focus mainly on ^Li production 
in these primordial QNe and extend the work of Ouyed (2013) by 
considering a more complete element composition in the (Pop III) 
type-II SN ejecta. We show that a QN can produce an ^Li plateau 
with values matching Spite levels of A(Li) ~ 2.2 and, with features 
satisfying the criteria set out above. 

The paper is organized as follows. In Section 2 we briefly 
present the QN model, primordial QNe and discuss their progen- 
itors. We give a brief overview of the QN model and of primordial 
QNe occurring in the wake of Pop III stars; i.e. QNe from the ex- 
plosions of massive neutrons stars left behind by core-collapse of 
stars at the lower mass end (2OM0 < Mprog. < 4OM0) of the Pop 
III mass distribution. Section 3 describes our model's parameters 
and its general features including ^''Ni destruction by QN neutrons 
and the resulting Fe impoverishment in the pristine clouds swept 
up by primordial QNe. Here we discuss the resulting sub-Fe spal- 
lation products and the a-elements and compare our findings to 
observations. When necessary the interested reader will be referred 
to Ouyed et al. (2011) and Ouyed (2013) for more details on the 
spallation process, numerical procedure, and relevant references. 
Section 4 focuses entirely on ^Li production and describes how a 
low-scatter, A{Li) ~ 2.2, plateau is obtained in our model. A discus- 
sion and predictions are given in Section 5 before we conclude in 
Section 6. 
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2 PRIMORDIAL QUARK-NOVAE 

The basic picture of the QN is that a massive neutron star (NS; 
with Mns > I.6M0) converts explosively to a quark star. Such an 
explosion can happen if the massive NS, in its spin-down evolu- 
tion or via mass accretion (e.g. from fall-back material), reaches 
the quark deconfinement density in its core (Staff et al. 2006) and 
subsequently undergoes a phase transition to the conjectured more 
stable strange quark matter phase (Itoh 1970;Bodmer 1971; Witten 
1984; see also Terazawa 1979). This results in a strange quark mat- 
ter conversion front that propagates toward the surface of the NS 
with a detonation occurring before the surface is reached (Keranen 
et al. 2005; Niebergal et al. 2010). The outcome of the QN explo- 
sion - besides the formation of a quark star - is the ejection of the 
NS's outermost layers (a very neutron-rich ejecta with an average 



mass Ml 



QN 



10 Mq) at relativistic speeds (with a Lorentz fac- 



tor averaging Fqn ~ 10; Ouyed&Leahy 2009). The outer layers are 
ejected from an expanding thermal fireball (Vogt et al. 2004; Guyed 
et al. 2005) which allows for ejecta with kinetic energy, £!SS > 10^^ 
erg. Hereafter we assume that a typical neutron in the QN ejecta 
will be streaming out with energy £qn ^10 GeV; not to be con- 
fused with the total ejecta's kinetic energy E^. This means that the 

Mqn ~ 10" M0 QN ejecta will have a typical total kinetic energy 
of£^~ 2x10^2 erg. 



2.1 Dual-shock Quark Novae (dsQNe) 

When the QN occurs a few days to a few weeks following the pre- 
ceding Type II SN (SNII) explosion, a dual-shock QN (hereafter 
dsQNe) is created. The interaction between the relativistic, neutron- 
rich, QN ejecta (Keranen et al. 2005; Ouyed&Leahy 2009) and the 
SNII ejecta leads to spallation of the SN material. This time delay, 
f(jelay, plays a crucial role in our model since it defines the density 
of the expanding SN ejecta (and thus spallation efficiency) when 
it is hit by the Eq^ ~ 10 GeV QN neutrons (the primaries in our 



model; see 5 3.1 below). If the time fdelay between SN and QN ex- 
plosions is lengthy, the SN ejecta will have dissipated such that the 
QN essentially erupts in isolation. However, when fdelay is on the 
order of days to weeks the neutron-rich QN ejecta interacts with the 
preceding SN ejecta leading to unique spallation products (Guyed 
etal. 2011; Guyed 2013). 



2.2 The Pop III Type-II SN (PopIII-SNII) progenitor 

Studies of star formation in the early Universe suggest that the 
first stars (Pop III) were typically more massive than present day 
objects, with lOOM© for zero metallicity (Z) objects (Abel et al. 
2002; Nakamura & Umemura 2002; Bromm&Larson 2004; Loeb 
et al. 2008). The progenitors of QNe are stars with masses in the 
2OM0 < Mprog. < 4OM0 range (e.g. Guyed et al. 2009a; Guyed et al. 
2009b; Leahy&Guyed 2009; Ouyed&Leahy 2012; Guyed 2013). 
Thus, primordial QNe are expected to occur in the wake of these 
Pop III stars from progenitors with mass range at the lower end of 
the mass distribution; 2OM0 < Mpiog. < 4OM0. These are expected 
to explode as Typell SNe and would most likely form NSs mas- 
sive enough to result in QNe once quark deconfinement density is 
reached in their cores (Staff et al. 2006). Hereafter we refer to them 
as PopIII-SNII progenitors. 



3 OUR MODEL 

In this section, we summarize and discuss our model's parameters 
and its general features and present some results. In all of our cal- 
culations we divide the target layer into A^^fpAr sublayers of radial 
thickness equivalent to a spallation mean-free-path (mfp), /isp.,AT- 
In each mfp of interaction, a given Aj nucleus will be hit multiple 
times by neutrons. To produce a realistic distribution of product nu- 
clei from this sub-layer, we draw the number of hits from a poisson 
distribution. As shown in Guyed et al. (201 1), a Poissonian descrip- 
tion is appropriate since we have a small spallation cross section but 
a rapidly increasing number of hits due to the cascading effect as 
more neutrons and protons are created. The interested reader is re- 
ferred to Guyed et al. (201 1; see eqs 3, 4 and 5 in that paper) and 
Guyed (2013) for more details on the spallation process, numerical 
procedure, and relevant references. 



3.1 Model Parameters 

We start by summarizing the key parameters in our model: 

• The PopIII-SNII ejecta's composition: Studies of the chem- 
ical composition of PopIII-SNII ejecta in the 2OM0 < Mprog. < 
4OM0 mass range can be found in the literature for mixed and 
unmixed SN ejecta (Umeda&Nomoto 2002; Nozawa et al. 2003; 
Cherchneff&Dwek 2009; Heger& Woosley 2010). Among the 
main products of these PopIII-SNII supernovae are ^^Ni,-'^S, 
2'*Si,2'*Mg,2''Ne,"'G and, ^^C. Other elements (including '^N; see 
03] below) are present at much smaller levels so we ignore them 
here. We also ignore elements heavier than ^^Ni in the present 
study. While some differences exist in mass yields between dif- 
ferent groups, and between models with mixing and non-mixing, 
we find that a typical 2OM0 PopIII-SN ejecta (with 10^' erg in 
energy) can be represented with mass yields Mnisn ~ 0. IM0, 
Ms,sN ~ O.O5M0, Msi.sN ~ O.O5M0, Mm^sn ~ O.O5M0, Mnc.sn ~ 
O.O5M0,Mo,SN ~ l-5M0,Mc,SN ~ O.15M0. We adopt these values 
as initial, fiducial, values for the target material; i.e. the PopIII-SN 
ejecta in our model. 

Although some mixing is likely to occur during the PopIII-SNII 
explosions, here we assume an onion-like profile of the expanding 
SN ejecta (i.e., no mixing) with the innermost ejecta, viz., ^^Ni nu- 
clei (mass number At=56) constituting the target Aj at a distance 
^At ('delay )=Vsn fdelay from the center of the explosion; we assume 
Vsn to be constant since the QN occurs while the SN ejecta is still 
in the Sedov phase (e.g., McCray 1985). Carbon ( C) makes up 
the last layer in this onion-like geometry overlaying the Oxygen 
layer. The '^He and H layers overlaying the carbon layer can be 
ignored since in most cases considered here we find that the spal- 
lation mean-free-path of the (primary, secondary and subsequent 
generations of) neutrons always exceeds the (H and ^He) layers' 
thickness. 

• The time delay (fdelay)^ The number density in the 
target layer (made of element Aj) is defined as wat = 
A/Ay,SN/(4w-Ki ARat^t^^h) where Mat.,sn is the mass of the layer 
(the target's initial mass in the PopIII-SN ejecta prior to impact by 
the QN neutrons), ARat the layer's thickness and /?At - I'sn'delay 
the radius of the layer from the center of the SN explosion; oth is 
the hydrogen mass. For a given PopIII-SNII explosion with known 
Ma^.sn (e-g- for 2 fixed amount of ^^Ni processed during the explo- 
sion) and given Vjn which remains constant early in the expansion 
of the SN ejecta, fdelay becomes the key parameter that defines the 
target density; i.e. the longer the time delay, fdelay the lower the 
target density when it is hit by the relativistic QN neutrons. 
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Table 1. Average multiplicity {() and energy of spallated neutrons for dsQNe with (delay = fouter ~ 8.3 days; this is the regime where the QN neutrons interact 
with the innermost layers (^*Ni, ^"S, ~*Si). For this time delay, A'mfpAT ~ 1 which means that the thickness of each target layer is equivalent to roughly one 
spallation mfp (see eq.ll|. Spallation stops at the -'^Si layer since the subsequent generation of neutrons exit the ^^^S layer with energies - 0.15 GeV which is 
below critical for spallation to ensue in the overlaying layers ^''Mg layer (Ssp-.Ms - 0.18 GeV; see discussion following eq. 2[. 



Primaiies (Eq = £qn = 10 GeV) Primaries (£0 = £qn = 5 GeV) 



At 



"Ni 



»Si 



"Ni 



»Si 



Sav.,i 


- 13.1 

- 0.77 
..44^; 

- 13.1 


-3.6 
-0.21 
-28Si 
-47.2 


- 1.4 
-0.15 
-2'Al 

-66 


- 11.3 
-0.44 
-«Sc 

- 11.3 


-2.8 
-0.16 

-28Si 
-31.7 


- 1.1 
-0.15 

-2'Al 
-34.8 


^ Apeak,! -M-iav.,i 
^net,i — ffav.,i-l ><^av.,i 



^The resulting spallation products would acquire a normal distribution peaking at atomic weight Apenk.i - Ax - ^,. 
* The net multiplicity (fnet; see eq. 4 in Ouyed 2013) translates to a total spallation neutrons and protons which amounts to M„+p = fnct^QN ~ 0.05-0.1 Mq. 



Table 2. Average multiplicity (f ) and energy of spallated neutrons in the oxygen ( '*0) layer for dsQNe with (delay > 'outer (here fdelay = 12 days); this is regime 
where the QN neutrons barely interact with the inner layers and proceed directly to the oxygen layer. The thickness of the oxygen layer (with Mq.sn = 1.5Mo) 
is equivalent to N^^ o - 10 sub-layers each of radial thickness or the order of a spallation mfp (see eq.FTj. Spallation in the oxygen layer ceases when the 
neutron energy is below the threshold value for O-spallation (iJsp.o - 0.27 GeV) or when the neutrons run out of target material (i.e. after Wsp.o - 10 mfps). 







Primaries (£0 = £qn = 


= 10 GeV) 




Primaries (£0 


= Eqn = 


= 5 GeV) 






mfp 1 


mfp 2 


mfp 3 


mfp 4 


mfp 5 


mfp 1 


mfp 2 


mfp 3 


mfp 4 


mfp 5 


fav..i 
p- £av..i-l (r-L^^J\ 


-3.8 
-2.7 

-3.8 


-2.7 

-1.0 

13c 

- 10.3 


-2 
-0.5 

14n 

-20.6 


-1.5 
-0.34 

14n 
-31 


-1.1 
-0.28 

-34 


-3.2 
- 1.6 

'3C 

-3.2 


-2.4 
-0.7 

14n 
-7.7 


-1.7 
-0.4 
14n 
-13 


-1.3 
-0.3 

-17 


- 1.1 
- 0.28 

- 19 





^The resulting spallation products would acquire a normal distribution peaking at atomic weight Apeak,] ~ 
* The net multiplicity (f„et; see eq. 4 in Ouyed 2013) translates to a total spallation neutrons and protons which amounts to Mn+p 



16 -f, 

= fnet^QN ~ 0.02- 0.04 Mq. 



The neutron spallation mean-free path (mfp) on a target Ax is 
/isp At = 1/(«At'''sp.,At) where the spallation cross section in milli- 
barns is (Tsp.^Aj = 45A*P mb (see Ouyed et al. 2011). The average 
number of collisions (i.e. the number of spallation mfps) an incom- 
ing neutron experiences in the target layer Aj is 

A^At „_ Mat,sn/0.1Mo 



A', 



mfp, At 



^0.69- 



'*sp.,At 



\03 



(^7-,56) (l'sn,5000fdelay,10)^ 



(1) 



where At,56 is the target atomic mass in units of 56 and 0. 1 Mq rep- 
resentative of the mass in the ^ Ni layer. In the equation above, I'sn 
is given in units of 5000 km s"' and is kept constant in all of the 
simulations shown in this paper. The time delay, /delay > is in units of 
10 days. For our fiducial values in mass of the target material and 
for ?delay = 10 days we get A'infp,Ni ~ 0.69,A'nifp,s ~ 0.41,A'nifp,si ~ 
0.42,A',rfp,Mg ~ 0.44,A'„fp,Ne ~ 0.47,iV„fp,o ~ 15.10, Afmfp,c ~ 1-64 
for spallation on ^''Ni, ^^s, ^^^i, ^'^Mg, 20]S[e, ISq, '^c, respec- 
tively. 

In a regime where /isp.,AT > ^^At (i-S- A'mfp,AT < 1)' there will be 
minimal or no spallation induced by the QN neutrons on the target 
Ax and the neutrons proceed onto the next, overlaying, layer. Thus, 
for a given Ma^.sn and Vsn this defines a critical time delay above 
which no spallation occurs in layer of element Aj. This allows us 
to define four regimes in our model: (i) /delay > ^outer with mini- 
mal or no spallation occurring in the inner layers. It means that 
for fdelay > 'outer ^ 8.3 days X (Mni.sn/O- lM0)l/2/(^a^5 ^,^^^^1^^^ 
days (see eq. 2 in Ouyed 2013), the QN neutrons barely interact 
with the inner layers and proceed directly to the outer (O and C 
layers); (ii) 'delay < 'inner with spallation occurring mainly in the 
inner layers, beneath the oxygen layer, with 'inner - 6.3 days x 



(MNi,SN/O.lM0)'/2/(AO-^5vs„,5OOo) (see Eq. 6 in Ouyed 2013^ 
This regime corresponds to the situation where the spallated (i.e. 
subsequent generations of) neutrons exiting the inner layers, do 
not have enough energy to spallate the outer CO-rich layers; (iii) 
'inner < 'delay < 'outer with spallation occurring both in the inner 
and outer layers. In this regime, the neutrons spallated in the in- 
ner layers end up with energies above the critical values for spal- 
lation in the '^O (and in certain cases in the '^C) to ensue; (iv) 
'delay > 'no-spall which effectively defines a regime of a QN oc- 
curing in isolation since no interaction occurs between the SN 
and QN ejecta. For our fiducial values, 'no-spall ~23-24 days. For 
'outer < 'delay < 'no-spalh Spallation would occur mainly in the oxy- 
gen and carbon layers. 

• The energy per nucleon, £qn, in the neutron-rich QN 
ejecta: It defines and sets the total multiplicity, f , of the spallated 
(i.e. subsequent generations of) neutrons and protons and their en- 
ergies. The average neutrons -l-proton multiplicity ^^y. on target Ax 
is given in eq. (2) in Ouyed (2013) as 



fav.-7AAT,56X(l+0.381n£„), 



(2) 



where £„ is the neutron energy in GeV which leads to a subse- 
quent generation of spallated neutrons and protons with energy ~ 
EnKm.- The resulting spallation product's atomic weight will peak 
at Ax - fav.- The <^av. > 1 condition gives the minimum neutron en- 
ergy necessary for spallation to occur in the Ni, S, Si, Mg, Ne, O, C 
layer to be E^^,^^-, ~ 0.105,0.139,0.153,0.173,0.206,0.268,0.416 
GeV, respectively. 

■* In Ouyed (2013), Pouter and finnei were defined as /ni and to, respectively. 
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Table [T] shows results from simulations involving the interaction 
of £qn = 5 GeV and Sqn = 10 GeV neutrons with the PopIII-SNII 
ejecta when fdday = 'outer - 8.3 days. For our fiducial values of the 
targets' initial masses this time delay implies Wnifp.Ni ~ Mnfp.S ~ 
^mfp.Si ~ ^mfp.Mg ~ ^mfp.Ne ~ 1 which means that on average a QN 
neutron will interact once per layer. An Eq = isQN =10 GeV QN 
(primary) neutron would lead to a multiplicity of ^av. ~ 13 in the 
^ Ni layer. This means a spallated (secondary generation of) neu- 
trons and protons exiting the layer with energy Ei = EqI^^m. ~ 0.77 
GeV. These secondary nucleons induce an average multiplicity 
^av. ~ 3.6 in the S layer which yields spallated (tertiary generation 
of) neutrons and protons exiting the layer with energy ~ 0.21 GeV. 
As can be seen from TablefT] spallation ceases at the Si layer since 
within the first three layers the neutrons energy has cascaded down 
to ~ 0. 15 GeV, below the critical value for Mg-spallation. The main 
(i.e. peak) spallation products are shown as well as the net multi- 
plicity fnet which yields a total mass of spallated neutrons and pro- 
tons which amounts to Mn+p = fnet^QN ~ 0.05-0. IMq. This table 
shows the pathway of a typical neutron with the expected average 
multiplicity given by eq. ([2|. However, because the f distribution 
is not uniform, some neutrons induce above average multiplicity 
(^ > ^av.) and stop within the ^^Ni layer while those inducing be- 
low average multiplicity (( < fav.) make it past the inner layers and 
induce spallation in the oxygen layer (see discussion below). 

Another example is shown in Table [2] which shows results from 
the interaction of £qn = 5 GeV and iiQN = 10 GeV neutrons with 
the PopIII-SNII ejecta when 'delay > 'outer (here 'delay = 12 days). 
For our fiducial values, A',^^ ,o ~ 10 and E^p.^o ~ 0.27 GeV. As can 
be seen from Table [2] neutrons interact with five sub-layers (each 
representing a spallation mfp) before their energy cascade down be- 
low 0.268 GeV. Thus half the amount of oxygen in the PopIII-SNII 
will be destroyed forming sub-0 products. In each sub-layer de- 
noted by index ; (i.e. for each mfp, or /lsp.,o)' 'hs atomic weight of 
the spallation products is close to a normal distribution peaking at 
atomic weight Apeak,! ~^T-^av.,i- The final abundances of the spal- 
lation products is a combination of these distributions (see Figure 1 
in Ouyed et al. 201 1). As can be seen in Table 2 the total spallation 
neutrons and protons amounts to Mn+p = (Tnet^QN ~ 0.02-0.04 Mq. 
This Table further illustrates the important point that most of the 
spallated neutrons exiting the oxygen layer will not be able to in- 
duce spallation in the overlaying C layer since their energy is below 
the threshold for C-spallation (~ 0.42 GeV). See however discus- 
sion below. 

We end this subsection by mentioning some important points: 

(i) The true ^ distribution is not uniform and is rather close to 
a normal distribution peaking at <fav.- This means that the spallated 
neutrons would exit a sub-layer ; of the target material (At) with 
a non-uniform energy distribution peaking at ~ £i/^av.,i- For the 
case shown in Table[2] most of the neutrons spallated in the oxygen 
layer will not be able to induce spallation in the overlaying C layer 
- since their energy is below the threshold for C-spallation £'sp.,C ~ 
0.42 GeV. However, neutrons leading to a multiplicity ^ < fav. do 
induc e spal lation after exiting the oxygen layer. As we show later 
(see 5 4.1.1 ', this has important consequences to production of ^Li 
from spallation in the C-layer; 

(ii) The logarithmic energy dependence of ^av. means that there 
is little difference between the multiplicity distribution induced by 
the Zsqn = 10 GeV and the iiQN = 5 GeV neutrons. However, for 
the lower energy primaries with Eq = Eq^ (i.e. QN neutrons), the 
resulting average energy per spallated nucleon is lower. In the case 
of spallation in the Oxygen layer for example (as illustrated in Ta- 
ble [2]l, the iiQN = 5 GeV simulations shows fewer neutrons with 



energy above £sp.,c ~ 0.42 GeV making it past the oxygen layer; 
i.e. spallation in the overlaying C-layer is reduced as compared to 
the iiQN = 10 GeV simulations (discussed further in 94. 1. U; 

(iii) Even if the spallation condition, E„ > £sp.,At' Is satisfied in 
any layer of target material Ax (a necessary condition), it is not suf- 
ficient when /isp.,AT > '^'^At^ i-s. when the spallation mfp exceeds 
the thickness of the target layer (see discussion following eq. [TJ. 
Thus the En > £sp.,At ^fl the A/?At > '^sp.Ax must be satisfied si- 
multaneously. 



3.2 - °Ni destruction and Iron-impoverishment in our model 

3.2.1 The t^g\.^y-[Fe/H] connection 

Following spallation in the ^^Ni layer, the resulting iron abundance 
([Fe/H]) in the primordial cloud swept by the combined SN+QN 
ejecta is (eq. 5 in Ouyed 2013): 



, 56 , Mnlsn/O-IMq 

- log nvi + log ■ 



■3.12. 



(3) 



where 77/ = M\/Mj^^sN is the normalized mass yields of prod- 
uct/element A following spallation on target Aj (originally present 
in the SN ejecta; here 56 refers to ^ Ni). Here Msw is the mass 
swept up by the mixed SN-l-QN ejecta in the primordial pristine 
cloud surrounding the site of the dsQN explosion; Ms,, is in units 
of 10^ Mq (e.g. Shigeyamac&Tsujimoto 1998; Machida et al. 2005). 

In our model, a dsQN with a given 'delay translates to a corre- 
sponding [Fe/H] abundance in the swept up cloud; this is set by the 
parameter 77^^ which is the level of depletion of the original PopIII- 
SNII ^ Ni following spallation. Figure 1 in Ouyed (2013) shows the 
heavy depletion of the original PopIII-SNII ^^Ni by spallation ob- 
tained for very short time delays ('delay < 2 days) where 7/5^ < 10" 
(i.e. leading to [Fe/H]< -7) . Hyper metal-poor stars (HMPs with 
[Fe/H]< -5 in the terminology of Beers&Christlieb 2005) would 
form from clouds swept up by 'delay < 3 days dsQNe in our model 
(see |5|. 

In dsQNe with long enough delays (i.e. 'delay > 'outer), there is 
minimal spallation in the Ni layer (i.e. rjt^ ^1) and in the inner- 
most layers ( S, S, Mg) layers. For such dsQNe, high values 
of [Fe/H] are obtained by increasing Mn[ sn or by decreasing the 
amount of swept up material, Msw (see eq.lSlabove). For example, 
a dsQN with /delay > 10 days, Mnjsn = 0.5Mq and M^w = lO'^'^Mg 
gives [Fe/H]~ -1.9. 



3.2.2 Suh-Fe elements 

Shown in FigurefTlin this paper are abundances of sub-Fe spallation 
products resulting from dsQNe with different 'delay ; shown is the 
relative abundance [X/Fe]. The original composition of the PopIII- 
SNII ejecta is shown by the plus signs in all of the panels. These 
elements are identified by an asterix next to the their names in all 
of the panels as their abundances vary with 'delay 

For 'delay > 'no-spall ~ 23-24 days (not shown here) no spalla- 
tion occurs and the original elements (identified with an asterix) fall 
exactly on the plus signs. For fouter ~ 8.3 days < 'delay < 'no-spalh 
there is minimal interaction in the inner layers and most spalla- 
tion products are light elements from O- and C-spallation. The 
original relative abundance of the PopIII-SNII elements is almost 
conserved except for some spallation in the oxygen layer which 
leads to the formation of B, Be and Li. For intermediate 'delay (1-6. 
'inner ~ 6.3 days < 'delay < 'outer). Spallation occurs in almost all of 
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the layers. As ?delay gsts shorter spallation efficiency start to de- 
crease in the outer layers and eventually for fdelay <~ 6 days spal- 
lation starts to become more prominent in the innermost layers. 
For even shorter fdelay (top panels), most of the QN neutrons are 
used up in the Ni layer thus shielding all of the overlaying layers. 
Effectively, the original PopIII-SNII sub-Ni elements (S, Si, Mg, 
Ne, O, C) present in the PopIII-SNII ejecta would be preserved in 
dsQNe with very short delays. In addition these elements will be 
produced as by-products of spallation in the ^^Ni layer which fur- 
ther increases their overall [X/Fe] values for fdelay < 3 days. The 
weak-dependency of multiplicity on the neutrons' energy explains 
the nearly similar final abundances for the £qn = 5 GeV and the 
iiQN = 10 GeV simulations. 



3.2.3 a-elements 

The so-called a-elements (O, Ne, Mg, Si, S, Ca, Ti) in very 
metal-poor stars are reported to be overabundant relative to iron 
at low metaUicity (e.g. McWiUiam 1997; Aoki et al. 2007) 
with [ajFe] > [tt/FeJo where [a/Fe] = ^([Mg/Fe] + [Si/Fe] + 
[Ca/Fe] + [Til Fe]) (the average abundance of Mg, Si, Ca and Ti). 
These enhancements increase linearly with decreasing metaUicity. 

Figure[2]shows the relative abundances of the original PopIII- 
SNII elements with respect to fdelay (top panel) and metaUicity 
(lower panel) for simulations with Fqn = 5 GeV. Figurelslshows the 
results for the Fqn = 10 GeV simulations. For fdelay > fouter ~ 8.3 
days, [X/Fe] remain close to the original values except for oxy- 
gen which experiences some spallation and thus depletion up to 
fdelay ~ 15 days since it is by far the most abundant PopIII-SNII el- 
ement. For finnei- < fdelay < fouter, spallation occurs in most of the 
PopIII-SNII layers leading to an overall reduction in the abun- 
dances of the original elements. This leads to a dip in the [X/Fe] 
values in Figure[2]and in FigurelSlfor -4 < [Fe/H] < -3. A plateau- 
like behaviour is noticeable in the dip for all of the PopIII-SNII 
elements. The [Si/Fe] and [Mg/Fe] levels in the dip region, com- 
bined with the copious production of Ti and Ca from Ni-spallation 
(see FigurefTJ lead to [a/Fe] > [Qr/Fe]o in our model and in partic- 
ular for fdelay < fouter ~ 8.3 days (i.e. [Fe/H] < -3) since spallation 
is more pronounced in the Ni layer than in the overlaying layers. 

An increase of [X/Fe] is evident for most elements when the 

metaUicity decreases below 4; i.e. when fdelay decreases below 

roughly 5 to 6 days. These layers are shielded at the expense of ^ Ni 
depletion explaining the increase in [X/Fe] much above the initial 
values. In general, and for short fdelay in our model, the more abun- 
dant a specific element is in the PopIII-SNII shocked ejecta, prior 
to spallation by the QN neutrons, the more it will appear abundant 
following the destruction of ^^Ni (i.e. the reduction of the final iron 
abundance in the swept up cloud). Since O and C are the two-most 
abundant elements in the PopIII-SNII ejecta, it is no surprise that 
[C/Fe] and [0/Fe] are very enhanced for very short delays (see be- 
low for the special of case of S). These enhancement would be re- 
flected in the composition of the gas cloud swept up by the SN-l-QN 
ejecta out of which the Pop II low-mass stars would form. 



3.2.4 Sulfur 

Figures [2] and [3] show that Sulfur is more depleted (with [S/Fe] 
showing the lowest plateau in the dip) than the other original ele- 
ments in the PopIII-SNII ejecta for 6.3 days < fdelay < 8.3 days (i.e. 
-4 < [Fe/H] < -3) due to its close proximity to the Nickel layer. 



Even for very short delays, because of the non-uniform distribu- 
tion in multiplicity <f, some spallated neutrons in the Ni manage 
to make it to the overlaying S layer This explains why the S 

plateau in the dip extends to metaUicity as low as [Fe/H] 4. 

Only for fdelay < 5 days (i.e. [Fe/H]< -4) does it become completely 
shielded and starts to follow the behavior of the other elements 
by increasing with decreasing [Fe/H]. ^^S has another distinctive 
feature: [S/Fe] exceeds [C/Fe] at very low metaUicity despite the 
higher initial abundance of C as compared to S in the PopIII-SNII 
ejecta. It turns out that Sulfur is among the main spallation products 
from Ni-spallation for fdelay < 5 days (i.e. [Fe/H] 4). This pro- 
duced S adds to the original S which explains why [S/Fe] exceeds 
[C/Fe] for our fiducial values. The trend of [S/Fe] versus [Fe/H] in 
extremely metal-poor stars requires further monitoring and obser- 
vations before its dependence on metaUicity can be confirmed (e.g. 
Nissen et al. 2007; Spite et al. 201 1 and references therein) and put 
our findings to test. 



3.2.5 Comparing to the average star 

Figure H] is a comparison of yields from a typical dsQN experi- 
encing spallation in the inner and outer layers (i.e. in the fi„ner < 
fdelay < fouter regime which corresponds to -4 < [Fe/H] < -3 in our 
model); we chose fdelay = 6.5 days. Our simulation results are com- 
pared to abundance pattern shared by the majority of the stars in 
the metaUicity range -4.0 < [Fe/H] < -3.0 with a normal [C/Fe] 
ratio (Cayrel et al. 2004; Spite et al. 2005; see also Figure 1 in 
Limongi&Chieffi 2012). As can be seen from Figure H] both the 
Fqn = 10 GeV and Fqn = 5 GeV simulations are successful in re- 
producing the observed abundances. Note the the high abundances 
of Fluorine and Scandium which are spallated in the innermost lay- 
ers in our model; these two rare elements are natural products of 
spallation (see §6.2 in Ouyed (2013). The Fqn = 5 GeV simula- 
tions underestimate the [C/Fe] and [0/Fe] abundances but this is 
easily remedied by increasing the original carbon and oxygen con- 
tent in the PopIII-SNII ejecta to reasonable amounts. 



3.3 Nitrogen in our model 

As shown in Ouyed (2013; see §3.3 in that paper), N forms by 
spallation in the A > 14 layers in the PopIII-SNII ejecta. How- 
ever, we find that '^N can also be enhanced by the injection of 
protons from spallation into the C-rich and 0-rich outer layers if 
temperatures in the PopIII-SN ejecta exceeds ~ 1.6 x 10^ K when 
it is hit by the QN ejecta. The carbon (and perhaps oxygen) could 
subsequently be processed into N via the CN(0) cycle (Amett 
1996). For an adiabatic expansion of the preceding SN shell (i.e. 
Fenv. oc f"^), the Fenv. > 1.6x 10^ K necessary condition for CN 
(proton capture) process to take place yields 



i?, 



fdelay <~ 5 days 



t1/2 
prog.,100Jg,„jo 

Vsn,5000 



(4) 



where i?prog.,ioo is the initial radius of the progenitor in units of 
100/?o (Rq is the solar radius) and Fenv.,io the initial temperature 
of the PopIII-SNII envelope in units of lO'" K. For our fiducial 
parameters, the above translates to (see Figure I in Ouyed 2013) 



(5) 



We find that the CN(0) processing of a small percentage (< 10%) 
of the carbon and oxygen into N is enough to account for the 
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Figure 1. Relative abundances [X/Fe] for spallation (sub-Fe) products versus mass number A (7 ^ A < 56; only stable isotopes are shown here and in other 
figures in this paper) for a dsQN with 5 GeV (left panels) and 10 GeV (right panels) QN neutrons. Letters show abundances from our simulations with the 
CN enhanced nitrogen (for short (jelay) shown as a filled diamond. The panels from top to bottom are for /delay(days) = 3,6,9, 15, respectively. Also shown are 
the resulting [Fe/H] and 'Li abundance, A(Li). The plus signs represents the initial relative abundance of the original PopIII-SNII elements; these elements 
are identified by an asterix (next to their names) as they evolve with ^deky Their initial abundances in mass are Mfji sn = 0. IM©, MssN = ^si.SN = ^Me.SN = 
^Ne.SN = O.OSMq, Mqsn = I.5M0 and Mc,sN = 0.15Mo. These initial abundances are the same in all of the figures in this paper unless stated otherwise. 
Similarly, unless stated otherwise, the mass of the cloud swept by a dsQN is set to Mjw = lO^M© in all of our simulations. 
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Figure 2. Relative abundances [X/Fe] versus ?deiay (top panel) and versus the con'esponding [Fe/H] (bottom panel) for ^^^S, -^Si, ^''Mg, ^"Ne, '*0 and '^C with 
initial relative abundances [X/Fe]x,SN=0.32, -0.01, -0.03, -0.28, 0.53, -0.09, respectively. These elements represent the original composition of the PopIII-SNII 
ejecta (i.e. the target material) in our model. Here, Eqn = 5 GeV. The horizontal dashed line in the bottom panel corresponds to [X/Fe]=0. 
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Figure 4. Relative abundances [X/Fe] of sub-Fe spallation products (elements identified by their names) versus mass number A from simulations with 
fdelay = 6.5 days; this resulted in [Fe/H] = -3.6. The CN enhanced nitrogen is shown by the filled diamond. Measured abundances of an average halo star 
(Cayrel et al. 2004; Spite et al. 2005) are also shown as plus {+) signs, with corresponding uncertainties, for a comparison. Left-hand and right-hand panels 
are for QN neutrons with iiQN = 5 GeV and 10 GeV, respectively. 



observed [N/Fe] values in an average star (in FigureHlthe CN en- 
hanced nitrogen is identified with a filled diamond). The C abun- 
dance should also be enhanced if CN equilibrium is established in 
which case the '^C/'-'C ratios should be of the order of a few (see 
however Caughlan & Fowler 1972). In our model, N production 
from CN cycling could occur at higher [Fe/H] when taking into 
account envelope reheating by the QN shock (following impact on 
the SN ejecta). However, we caution that our model is more detailed 
than the standard CN cycle and would need further investigations 
for better estimates of '^N enhancement and the expected '^C/^^C 
values. 



4 LITHIUM SYNTHESIS IN OUR IVIODEL 

The ^Li produced by spallation on a target Aj in the ejecta of the 
PopIII-SNII ejecta is given by eq.(14) in Ouyed (2013): 



A(Li) - logrif^^ + log ^ ° + 6.33 , 

' M,w/105Mo 



(6) 



„At 



where J/, ^ = Mii^JMp^^^^ is the normalized mass yields of Li 
following spallation on target Ax with initial mass /Wat,Sn; the 
1.5Mq is representative of our fiducial value for the oxygen con- 
tent in mass in the PopIII-SNII ejecta. The final ^Li abundance is a 
combination of ^Li produced from spallation on each of the targets 
present in the PopIII-SNII ejecta. 

Shown in Figureplis the ^Li abundances resulting from spalla- 
tion by iiQN = 5 GeV (left panels) and Zsqn = 10 GeV (right panels) 
QN neutrons on a typical PopIII-SNII ejecta in our model. The top 
panels show A(Li) from simulations spanning a range in ^^Ni con- 
tent of the PopIII-SII ejecta {O.OSMq < A/ni.sn ^ 0.5Mq) and a 
range in /delay (2 days ^ fdelay ^ 24 days). The bottom panels show 
the corresponding metallicity [Fe/H]; as explained in |3.2.1| there 
is a direct connection between fdelay and [Fe/H] in our model. For 
fdelay > 24 days the QN neutrons traverse the PopIII-SNII ejecta 
with minimal interaction; i.e. no spallation occurs and the PopIII- 
SNII ejecta abundances would not be modified or affected (i.e. 
T]J = I). While Ni content was varied, the other initial abun- 
dances in mass in the PopIII-SNII were set to their fiducial values: 

Ms,SN = Msi,SN = MMg.SN = MNe,SN = 0.05 /Wq, /Wo,SN = I-SMq 



and Mc^SN = 0.15Mq. We adopt a 10^ Mq of pristine cloud (the 
cradle of the Pop II low-mass stars) swept up by the SN-l-QN ejecta 
(e.g. Shigeyama&Tsujimoto 1998). 

The different panels in Figure[6](Figure[8]) show the resulting 
A(Li) versus /delay for each target for the Eqn = 5 GeV (£qn = 10 
GeV) simulations. The corresponding A{Li) versus [Fe/H] is shown 
in Figure [7] (Figure [9]l. The lower right panel in Figure [6] (Figure 
[Sl shows the combined A(Li) versus /delay which effectively corre- 
sponds to the projection of the upper left (right) panel in Figure [5] 
onto the A(LO-/delay plane. Similarly, the lower right panel in Fig- 
ure [7] (Figure [9) shows the combined A(Li) versus [Fe/H] which 
effectively corresponds to the projection of the lower left (right) 
panel in Figure [5l onto the A(Li)-[Fe/H] plane. One can see from 
these panels that^Li is produced in the ^^Ni and •'^S layers when 
'delay < 'outer and in the oxygen and carbon layers for longer delays 
(/outer ~ 8.3 days is shown by the vertical dashed lines in the fig- 
ures). No ^Li is spallated in the intermediate layers (Si, Mg, Ne) 
and in general ^Li is mainly produced either in the innermost lay- 
ers or in the outermost CO layers. The vertical doted lines defines 
the /Vuifpc = Ai?c/'lsp.,C = 1 boundary which defines a timescale 
above which the spallation mfp exceeds the thickness of the C- 
layer; i.e. no C-spallation occurs for longer time delays even if the 
neutrons exiting the O-layer have energies exceeding £sp.,C ~ 0.42 
GeV which is above the critical value for C-spallation to ensue. 

There is a gap in the ^Li abundance (i.e. a drop in the A(Li) 
value) in the intermediate regime, /i„ner < 'delay < 'Ni (which corre- 
sponds to -4 < [Fe/H] < -3) since most neutrons end up being 
used for spallation of A > 7 elements in the inner layers. Thus, 
for Pop II low-mass stars formed from clouds swept up by dsQNe 
with /delay < 8.3 days (i.e. for [Fe/H]< -3), the resulting ^Li abun- 
dance would show a larger scatter mainly because in these dsQNe 
^Li is not the main spallation product and the neutrons are used 
up to make many other spallation products. At this point, we must 
note that the 'Li abundance in dsQNe leading to [Fe/H]< -3 (i.e. in 
'delay <~ 10 days dsQNe) should be considered upper limits in our 
model. For such dsQNe, the SN ejecta is still hot {T > 2.5 x lO*" K) 
when it is spallated. Combined with the copious production of pro- 
tons (see Tablesll|andl2|(, we expect most of the spallated ^Li to be 
destroyed by reactions with the protons or at least become reduced 
(see §5 in Ouyed 2013). 



10 Rachid Ouyed 



4.1 The high-scatter, A{Li) ~ 2.2, plateau 

Two "plateaus" can be seen in the top right panel in Figure IS] (i.e. 
for the iiQN = 10 GeV simulation). One resulting from O-spallation 
with 1.8 < A{Li) < 2.5 and a second one from C-spallation with 
2.5 < A{Li) < 3.5. These two plateaus can also be seen in panels 
labeled "C-Target" and "0-Target" in Figures [S] |7|8[ [9] Both oc- 
cur for fdelay > fouter ~ 8.3 days which means that the QN neutrons 
proceed directly into the outer layers with minimal interaction with 
the Ni, S, Si, Mg and Ne layers. Spallation in the 0-layer leads to 
smaller ^Li abundance and more scatter in the ^Li plateau since the 
neutrons spallate more products than in the case of spallation in the 
C-layer. 

We first focus on the A{Li) ~ 2.2 plateau uncovered in our 
simulations (the A{Li) = 2.2 value is indicated by the horizontal 
planes in FigureBland by the horizontal lines in Figures [6] [7] [8] [9]l. 
This plateau is populated by ^Li spallated from oxygen for dsQNe 
with fdelay > 8-3 days. As can be seen from Figure |5| while the 
Eqn = 5 GeV simulations yield A(Li) approaching the 2.2 value 
from the lower bound, barely exceeding it, the Eqj^ = 10 GeV sim- 
ulations yield ^Li abundances oscillating around A(Li) ~ 2.2 (with 
1.8 < A(Li) < 2.5). To understand this interesting finding from our 
simulations we first note that in each sub-layer ; (i.e. for each mfp, 
'isp.,o) of the 0-layer, the distribution of the atomic weights of 
the spallation products is close to a normal distribution peaking at 
^peak.i ~ Aj - (.^^ { where fav.,i is the average multiplicity for sub- 
layer ( (see Table[2j. For each sub-layer the corresponding amount 
in mass of spallated ^Li can be estimated from the corresponding 
distribution. 

Since the highest ^av. is induced in the first sub-layer this 
leads to a peak spallation product (Apeak = IS-^Tav.) closest to ^Li. 
Subsequent sub-layers would yield peak spallation products get- 
ting closer to the target's Aj (here Aj = 16) thus reducing the ^Li 
contribution to the overall distribution of the spallated elements 
and thus a reduction in the final amount of ^Li spallated. How- 
ever, for subsequent sub-layers the multiplicity is higher and more 
neutrons+protons induce spallation. The final abundances of the 
spallation products is a combination of these distributions (see e.g. 
Figure 1 in Ouyed et al. 2011) which yields Mli ~ ffy^netA^QN - 
PjMqj^ where /3j = ffy^net- Since the ^Li contribution given by 
the distributions is roughly aj ~ 0.01 (of the order of ~ 1%) and 
(net ~ 10 (see e.g. Table[2]( this yields /?7 ~ 0.1 The corresponding 
^Li abundance A(Li) = logMLj/A'^ + 12 is then roughly estimated 
to be 



A(LO~ 2.15+ log L87,().l 



Mqn/10-^M, 
M.w/lO^Mo 



(7) 



where we made use of Ml; = Mii^Jlmn and Nu = Msvilmn. The 
amount of ^Li spallated is roughly equivalent to a normalized mass 
yield zy'^ ~ 6.7 x 10"^ or logrjl!' — 4. 17 for our fiducial values; i.e. 
A{Li) ~ 2.15 from using eq.jSl. 

The weak dependence of the multiplicity on energy leads to 
similar distributions of spallation products in the oxygen layer for 
the range in £qn expected in QNe explosions. This translates to 
similar ^net and a-; (i.e/S?) values which explains why the £qn = 5 
GeV also hovers around the same A(Li) value and effectively im- 
plies that A(Li) is mainly dependent on Mqn and M^v,. In fact by 
scaling Mj^ with Mqn (i.e. with the kinetic energy, E^, of the 
QN ejecta) we always recover an A{Li) close to Spite values when 
fdelay > foutei-; 1-6. when spallation occurs mainly in the outer oxygen 
and carbon layers as shown in ^4.2|below. 



4.1.1 The C-spallated Li plateau 

When Ai?c < 'isp.,C> sven if enough neutrons with energy above 
£sp.,C ~ 0.42 GeV make it past the O-layer, they will not inter- 
act with the overlaying carbon and there will be no spallated ^Li 
and the corresponding plateau would disappear. From eq. (fTll, the 
Mnfp.C > 1 condition (i.e. /S.Rc > /isp.,c) yields fjelay <~ 13 daysx 
(Mc,sn/0. ISMq)'^ for our fiducial values. This means that spal- 
lation in the C-layer would occur mainly when fouter ~ 8 days < 
'delay < 13 days which is a narrow window. On the other hand, if 
we increase the carbon content in the PopIII-SNII ejecta as to make 
A/?c > 'isp.,C then spallation would ensue. Figure 10 shows the two 
plateaus again for £qn = 10 GeV but for much higher carbon con- 
tent than our fiducial values (here Mc,sn = 'Wq.sn = 1.5Mq). As 
can be seen in the figure, although the plateau form O-spallation is 
close to the Spite value, it is clearly dwarfed by the ^Li produced in 
carbon layer. For Mc^sn = 1.5M0, the ARc < 'isp.,C regime occurs 
when fdelay > 30 days beyond which the C-spallated ^Li plateau 
disappears. We find that in PopIII-SNII with high carbon content, 
and in particular those with iiQN ^10 GeV, the ^Li plateau from 
C-spallation is dominant and far exceeds the Spite value. Another 
possibility would consist of reducing oxygen content in the PopIII- 
SNII ejecta so that the primary (i.e. QN) neutrons experience only 
a very few interactions (mfps) in the oxygen layer before making 
their way towards the C-layer. This is the least likely scenario since 
it seems that oxygen is produced in much larger quantity than car- 
bon in PopIII-SNII ejecta (e.g. Umeda&Nomoto 2002). 

On average, spallated neutrons in the Eqn = 5 GeV case ac- 
quire lower energies than the Eqn = 10 GeV case. Thus, there is 
no substantial amounts of ^Li spallated in the C-layer since fewer 
> 0.42 GeV neutrons manage to make past the O-layer. The result- 
ing plateau is very scattered and remains below Spite values; see 
the panels labeled "C-Target" in Figures|6]and|7] 



4.2 A low-scatter, A{Li) ~ 2.2, plateau 

If we associate high Eqn dsQNe with short time delays and low 
iiQN dsQNe with relatively longer delays (> 15 days), a low-scatter 
plateau results as shown in Figure [TT] The top panel shows abun- 
dances from simulations with Eq^ = 10 GeV and Ms„ = 10^ Mq 
with 5 days < ?delay < 15 days. The middle panel shows results 
from an £qn = 5 GeV simulation with Msw = 10 Mq and for 
15 days < fdelay < 30 days. The combined ^Li abundance from these 
simulations leads to a plateau with features reminiscent of the ob- 
served one. To extend the plateau to higher [Fe/H] values we had to 
include simulations where the mass of the pristine cloud swept up 
by the dsQN is lower than the fiducial value of 10 Mq. For exam- 
ple, the Msw = 10 Mq simulations effectively shift the plateau to 

the right in Figures[7]and[9]by yielding [Fe/H] values close to 2 

(see eq.lSl. However, a decrease in Msw increases the A{Li) value 
(see eq.pb which is remedied by a decrease in £qn which decreased 
the amounts of spallated ^Li accordingly. The high A(Li) ~ 3.5 val- 
ues in the top panel correspond to abundances from C-spallation. 
However, the thinnest plateaus result from simulations which min- 
imize spallation in the C-layer. This meant setting Mc,SN < O.IMq 
for our fiducial values, or in general by taking Mc,sn < A/q.Sn/IO. 
While more sophisticated calculations may be required to fur- 
ther constrain the parameters, the QN model naturally explains this 
dichotomy: In general, the reduced mass-loss in Pop III stars means 
that the more massive progenitors could in principle end up with 
the more massive cores at explosion and should thus lead to heav- 
ier NSs. Or, if fallback is greater from the more massive QN pro- 
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genitors (i.e. those closer to the black hole hmit which we set here 
to be 40Mo), this should also translate to the formation of a more 
massive compact remnant (i.e. a parent NS with higher core den- 
sity). These heavier NSs should undergo a QN explosion on shorter 
time delays as found by Staff et al. (2006). The core-collapse of the 
massive progenitors followed by the QN explosion of the massive 
NS would likely lead to a more energetic SN-l-QN ejecta which 
should sweep up more pristine material (i.e. Ms„ should be higher 
in these dsQNe; e.g. Shigeyama&Tsujimoto 1998). Naturally, the 
shorter time delay (i.e. the lower [Fe/H]) dsQNe would lead to 
a higher scatter in the plateau, and its eventual meltdown for the 
shortest delay dsQNe, since ^Li will be produced in the inner lay- 
ers from multiple targets. In this picture, low-mass Pop II stars 
in the Spite plateau would have originated from pristine clouds 
swept up by dsQNe with progenitors in the lower end of the mass 
range, 20Mq < Mpiog. < 40Mq, with reduced fallback. As discussed 
later (see ^, dsQNe with short delays should be the rarest among 
dsQNe. The implication is that the resulting Pop II stars (those with 
very low [Fe/H] and high A(Li) scatter according to our model) 
would be the rarest among very metal-poor stars. 



4.3 Discussion 

We now discuss some implications of our findings to Lithium in the 
Universe: 



• The BBN component: In this work and in Ouyed (2013), we 
have neglected the BBN contribution in the cloud swept up by the 
SNh-QN ejecta and considered only A(Li) from dsQN spallation on 
PopIII-SNII ejecta. Among plausible scenarii that could lead to the 
destruction of the BBN component in pristine clouds, before they 
get swept up by the combined SN+QN ejecta, are: 

(i) The QN provides copious amount of neutrons with some of 
them freely streaming ahead of the combined SN-fQN ejecta. These 
free QN neutrons could destroy the BBN ^Li (e.g., see Albomoz 
Vasquez et al. 2012 and references therein for relevant nuclear re- 
actions). This might contribute to the destruction of the surrounding 
BBN ^Li before material is swept up by the SN+QN ejecta - this 
remains to be confirmed; 

(ii) A massive astration of ^Li by processing in Pop III stars as 
suggested by Piau et al. (2006). This model provides a natural con- 
nection to ours since it appeals to Pop III stars in the IOM0-4OM0 
mass range; the favored QN progenitor mass range. Thus, adopting 
the Piau et al. (2006) hypothesis, one could then envision a sce- 
nario in which the BBN ^Li is processed and destroyed by the Pop 
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Figure 9. Same as in figurelTlbut for the Eqn = 10 GeV simulations. 
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III star before it undergoes a PoPIII-SNII explosion and the sub- 
sequent QN event. In this scenario, the issue of overproduction of 
CNO elements plaguing the Piau et al. (2006) hypothesis is alle- 
viated since in our scenario these would be reduced by spallation. 
The pristine material/cloud swept up would then carry the imprint 
of ^Li resulting from spallation of PopIII-SNII ejecta by dsQNe 
without the overproduction of CNO elements; 

(iii) If it turns out that neither of these scenarios are satisfactory, 
we might envisage the bold but very unlikely proposition that ^Li 
might not have been processed during BBN and is rather a relic of 
primordial dsQNe spallation. 

• Super-Li rich candidates: A few super-Li rich stars have 
been reported in the literature exceeding the BBN value (e.g. 
the A{Li) - 2.87 star in the globular cluster M4; Monaco et al. 
2012). Other stars far exceed the meteoritic values of A(LO = 3.25 
(Deliyannis et al. 2002; Koch et al. 2012; Adamow et al. 2012). 
Standard enhancement scenarios including planetary ingestion, Hy- 
pernovae, outward radiative acceleration of ^Li from deeper regions 
in the stars, and binary transfer face challenges in explaining such 
extreme levels of Li-enhancement (e.g. Monaco et al. 2012; Koch 
et al. 2012). 

In our model, such high A(Li) values result from spallation in 
the C layer. C-spallation would occur for a narrow window of time 
delay (8.3 days < fdelay < 13 days for our fiducial values) and only if 
enough > 0.42 GeV neutrons make it past the O-layer; i.e. mainly 
in dsQNe with £qn ^10 GeV and with Mc,sn > 0.1Mo,SN (see 
discussion in ^4.1[ (. These stringent constraints suggest that super- 
Li candidates (i.e. Pop II low-mass stars with ^Li spallated in the 
C-layer in our model) should be relatively rare compared to normal 
(i.e. with A{Li) ~ 2.2 from 0-spallation) candidates. 

• Interstellar ^Li abundance: Howk et al. (2012) reported 



measurements of interstellar ^Li abundance of the Small Magel- 
lanic Cloud with values nearly equal to the BBN predictions. This 
suggests that contributions of galactic chemical evolution can be 
neglected in this case. We suggest instead that the measured value 
may be a signature of a gas swept up by a dsQN with important 
C-spallation (see discussion above). Interestingly the Li/^Li=0. 13 
measured by Howk et al. (2012) is very close to the value predicted 
by Ouyed (2013; see §5.3 and Fig. 8 in that paper). The high Li/^Li 
values are expected in our model since ^Li and ^Li would be pro- 
duced with nearly equal proportions from spallation. 



5 DISCUSSION AND PREDICTIONS 

Below we discuss some implications of our model to the early uni- 
verse followed by some specific predictions: 

• Formation of Pop II low-mass stars: There are two models 
for the transition from Pop III to Pop II discussed in the litera- 
ture: (i) Atomic fine-structure line cooling (Bromm & Loeb 2003). 
In this model, low-mass star formation can occur in gas that is 
enriched beyond the critical abundances of [C/H]crit — 3.5 ± 0.1 
dex and [0/H]crit ~ -3 ± 0.2 dex (or Z > Z^. ~ IO^^'^Zq); (ii) 
Dust-induced fragmentation (e.g. Schneider et al. 2006a). The dust- 
cooling model predicts critical metallicity, Z, smaller by a factor of 
10-100 than those predicted in line cooling models. 

The discovery of SDSS J102915 + 172927 (Cafi'au et al. 2011), 
an ultra metal-poor (UMP) star with [Fe/H] 5 (and with no ex- 
cess of carbon) bom from a cloud with Z < Zcii challenges the line 
cooling model (see also Klessen et al. 2012) and seems to hint at a 
dust-induced fragmentation. However, it has been argued that dust 
production in very metal-poor gas is too inefficient (e.g. Silvia et al. 
2010; see however Nozawa et al. 2012) in which case one should 
consider fragmentation at high densities as another plausible mech- 
anism to form low-mass stars in the pristine universe (e.g. Greif et 
al. 2011 and references therein). In our model, Fe-poor does not 
translate to Z-poor since in short delay dsQNe the outer CO-rich 
layers should be preserved. The swept up cloud will be Fe-poor but 
not necessarily Z-poor. One can also have the situation where most 
of the PopIII-SNII C and O is spallated without destruction of the 
innermost layers. Furthermore, the extremely dense relativistic QN 
ejecta should provide high enough compression to drive the needed 
fragmentation in the pristine cloud swept-up by the dsQN. The 
cloud could in principle fragment into dense cores whose masses 
should be comparable to the Jeans masses and collapse to form low- 
mass stars that can survive to the present (Machida et al. (2005). 

Figure [12] compares relative abundances of stable isotopes in a 
pristine cloud swept up by a /delay = 3.5 days dsQN to measured 
values in SDSS J102915-I-172927. Both the £qn = 5 GeV and the 
iiQN = 10 GeV simulations give very encouraging fits to measured 
abundances although our model overproduces Ti. As noted ear- 
lier, large abundances of Fluorine and Scandium are expected in 
our model since these rare nuclides are a natural outcome of spal- 
lation processes (see §6.2 in Ouyed 2013). The ^Li shown are up- 
per limits in our model since as we have noted, for fdelay <~ 10 
days (i.e. for [Fe/H]< -3), the PopIII-SNII envelope temperature 
exceeds 2.5 x 10* K when spallation occurs which should reduce 
the final abundance of the spallated ^Li. 

• The primordial, Pop III, IMF: Our model requires no mod- 
ification to the primordial IMF (the very massive Pop III found in 
simulations; e.g. Wise&Abel 2007) and does not contradict obser- 
vations that hint at contamination of Pop II star-forming clouds by 
the explosions of 10Mq-50Mq progenitors (e.g. Umeda&Nomoto 
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Figure 11. Top panel: A(Li) versus [Fe/H] for simulations with £qn = 10 GeV, M^„ = 10^ Mq and for 5 days < (delay < 15 days. Middle panel: A(Li) versus 
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2005; Tumlinson et al. 2004). In a heavy Pop 111 IIVIF, stars with 
2OMO-4OM0 would lead to dsQNe. Thus our model might help 
reconcile the Pop 111 simulations and the observations of Pop 11 
low-mass stars if the latter formed from pristine clouds swept-up 
by primordial dsQNe as described in this woris (see §6.3 in Ouyed 
2013 for more discussion and relevant references). 

• The statistics; In a heavy Pop 111 IIVIF (e.g. see discussion in 
Larson 1998; see also Schneider et al. 2006b and references there 
in), dsQNe progenitors would contribute a small percentage of the 
population. Thus, according to our model, the statistics of the re- 
sulting low-mass Pop 11 stars in the galactic halo may be a reflec- 
tion of the IIVIF of Pop 111 stars (specifically the lower mass end 
of the distribution). Interestingly, the statistics of the Icnown old, 
[Fe/H] < -2 stars are intrinsically rare in the Galaxy and represent 
only a small percentage of the total stellar mass (see Table 3 in 
Beers&Chiisdieb 2005). 

Furthermore, in the QN model, the more massive the parent neu- 
tron star (to be experiencing a QN explosion), the more likely it will 
reach deconfinement in short time delays (Staff et al. 2006). As we 
have mentioned earlier, if fallback or the size of the collapsing core 
increases with the QN progenitor's mass, then these should lead 
to heavier NSs which are more likely to undergo a QN transition 
on shorter time delays. This means that short fdelay dsQNe should 



be the rarest in a heavy Pop 111 IIVIF peaking below ~ 40Mq. In 
other words, the lower the [Fe/H] (i.e. the shorter the fdelay) the 
smaller the statistics of the ensuing Pop 11 stars. In this picture. 
Hyper- IVIetal Poor (HIVIP) stars (with [Fe/H]< -5 in the terminol- 
ogy of Beers&Christlieb 2005) would have originated in pristine 
clouds swept up by dsQNe from progenitors on the high end in the 
20Mq < Mpios. < 40Mo mass range. However, this conclusion re- 
lies heavily on the exact shape of the Pop 111 DVIF which remains to 
be confirmed. 

• Origin of HIVIP stars; We focus here on two known HIVIP 
stars, HE0107-5240 ([Fe/H] = -5.54; Christlieb et al. 2002) and 
HE1327-2326 ([Fe/H] = -5.76; Frebel et al. 2008). These two 
stars have extremely high enhancements of the light elements C, 
N, and O relative to Fe. Figure [T3] show a comparison of the pecu- 
liar chemical composition of HE 0107-5240 and HE 1327-2326 
to abundances from dsQNe with for fdelay = 2.8 days and fdelay - 2.6 
days, respectively. The measured relative abundances of elements, 
[X/Fe], in these sources we take from Norris et al. (2013; see their 
Table 4). In our model, HlVIPs would correspond to the scenario 
where substantial interaction occurs in the inner layer either using 
up the majority of the QN neutrons or reducing their energy be- 
low the critical value for O- and C-spallation to ensue (see §3 in 
Ouyed 2013). The shielding (from spallation by the QN neutrons) 
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of the outer CO layers combined with the substantial destruction 
of the Ni layer explains the large [C/Fe] and [0/Fe] in these 
objects (see Figures [T] [2] and [Sjl. We find a good fit to measured 
abundances except for Ti which is overproduced in our model. 
There is a strong apparent depletion of lithium, (A{Li) < 1.5) in 
these two HMP stars. A strong lithium depletion in the atmosphere 
of HE1327-2326 is unlikely given its mass and temperature Teff ~ 
6000 K. The more massive HE0107-5240 may have depleted its 
own ^Li but it is hard to imagine a process that would practically 
bring its abundance much below the plateau. Our model leads to the 
right levels of ^Li abundances although the ^Li abundance from our 
simulations should be considered upper limits since [Fe/H]< -3 for 
HMPs (see pO] ). 

There have been investigations on the origins of the HMP stars. 
It has been argued that the binary mass transfer hypothesis provides 
best fits to the overall abundances of these stars (Suda et al. 2004). 
This hypothesis remains to be confirmed by the discovery of ra- 
dial velocity variations in HMP stars. Furthermore, this hypothesis 
implies that the stellar IMF at Z < IO^^Zq was biased towards in- 
termediate mass stars (Tumlinson 2007). The large C enhancement 
with respect to Fe in HMPs stars we argue, might not be an indi- 
cation of binarity (e.g. Suda et al. 2004) but rather a direct conse- 
quence of important Ni (and thus Fe) destruction and a preserva- 
tion of the carbon in the PopIII-SNII ejecta. In the binary accretion 
scenario, high level of C should be accompanied by high levels of 
s-process which is not always observed (e.g. Cohen et al. 2006). In 
the dsQN model, slow neutrons capture processes and the resulting 
s-process products become important in dsQNe with /delay leading 
to HMP stars (see §4 in Ouyed 2013 for more details; in particular 
eq. 1 1 in that paper). We note however that the extreme N abun- 
dance in HE1327-2326 can be better fit in our model if we assume 

N to be present in reasonable amounts in the PopIII-SNII ejecta. 
Interestingly, it has already been argued that '"^N could be pro- 
duced in Pop III stars close to the ~ 4OM0 limit (Heger&Woosley 
2010); this limit is lower when taking rotation into account (Jog- 
gerst et al. 2010). In our proposed picture, the extreme [C/Fe] and 
[0/Fe] values in HE1327-2326 is already indicative of a dsQN pro- 
genitor with mass at the high end of the 20Mq < Mprog. < 4OM0 
mass range. It is only fitting and natural to assume some ''^N al- 
ready present in the PopIII-SNII ejecta. Of course, our model does 
not rule out the possibility of a companion at play in the case of 
HE1327-2326 given the initial CNO abundances we had to adopt 
for the fits. 



5.1 Predictions 

We now list some predictions of our model: 

• An immediate implication of the astration-dsQN scenario sug- 
gested in this paper (see discussion in 5 4.3 1 is that the BBN ^Li 
should be preserved in environments surrounding > 4OM0 Pop 
III stars. Pop II low-mass stars formed from clouds swept up by 
> 40Mq PopIII-SNII should show imprints of elements synthe- 
sized in these massive Pop III stars and carry a BBN ^Li abundance 
(in particular in warm low-mass stars formed these clouds). These 
stars would show a BBN Li/^Li ratio of ~ 10" unlike the low- 
mass stars formed in clouds swept-up by dsQNe where, besides the 
A{Li) ~ 2.2 abundance of ^Li from 0-spallation, much higher val- 
ues of Li/^Li are expected (see §5.3 and Figure 8 in Ouyed 2013). 
In addition, objects born from pristine clouds swept up by dsQNe 
should bear chemical imprints of spallation products as described 



in this work. Thus we predict two distinct classes of low-mass Pop 
II stars. 

• We found that in dsQNe with 8.3 days < /delay < 13 days and 
£qn Js 10 GeV, spallation in the C-layer resulted in ^Li abundances 
exceeding meteoritic values {A{Li) ~ 3.25). The destruction of C at 
the expense of lighter elements in these dsQNe means that super-Li 
rich low-mass stars should also show relative depletion in Carbon. 
However, C is also a by-product of spallation in the O layer 
which suggest that not all super-Li candidates should appear carbon 
poor. 

• In our model, QNe with fdelay < 2 days would lead to exces- 
sive (almost complete) destruction of Ni. Thus, low-mass Pop II 
stars with much lower [Fe/H] values than HMP stars should even- 
tually be discovered if QNe with extremely short time delays occur 
in nature. However, these dsQNe are statistically extremely rare 
since they would originate from stars lying the closest to the black 
hole limit (~ 40Mo): those forming the most massive NSs which 
undergo a QN phase transition with the shortest time delays. 

• If some low-mass Pop II stars are made in binaries in the pris- 
tine cloud swept up by dsQNe with very short delays (/delay < 2.5 
days or [Fe/H]< -6), these should show excessive enhancement of 
CNO elements. The preservation of the original C and O (and N 
for the more massive dsQN progenitors) elements present in the 
PopIII-SNII ejecta, their production by spallation in the inner lay- 
ers, and their accretion from the companion should lead to [C/Fe], 
[0/Fe] and [N/Fe] values much in excess of those measured so far 
in carbon-enhanced, extremely metal-poor (CEMP) stars. The dis- 
covery of such stars, in particular if they show large *Li/^Li ratios 
and large abundances of specific spallation products such as Flu- 
orine and Scandium (see e.g. Figures [T] [12] and |13^ , would be a 
strong indication that our model might be a correct representation 
of extremely metal-poor Pop II stars and of the origin of the Spite 
plateau. 

• We echo a previous point that besides being very neutron-rich, 
the QN ejecta is also rich in A > 130 elements (Jaikumar et al. 2007; 
see discussion in §6.4 in Ouyed 2013). A typical QN ejects on av- 
erage ~ 10 -10" Mq of ejecta rich in A > 130 elements which is 
10-100 times larger than in the case of a typical core-collapse SN 
ejects where on average 10" Mq in heavy elements is ejected (e.g. 
Woosley et al. 1994 and references therein). Since we estimate that 
on average QNe occur at a rate of 1 for every 100 core-collapse SNe 
(Leahy&Ouyed 2009; Ouyed et al. 2009a&b), this means that con- 
tamination by QNe should be equally important if not more domi- 
nant. These elements should be detectable in environments around 
low-mass Pop II stars with ^Li abundance reflective of dsQNe. A 
possible site for the QN A > 130 tracers would be the ISM envi- 
ronment studied by Howk et al. (2012) in the Small Magellanic 
Cloud where the super-Li A{Li) value and the high Li/^Li mea- 
sured, as we argued previously, might be indicative of a dsQN with 
C-spallation imprint. 

• In this work and in Ouyed (2013), we argued that extremely 
metal-poor Pop II low-mass stars could have been bom in pristine 
clouds swept up by primordial dsQNe. Furthermore, in Ouyed et 
al. 2009a we showed that primordial dsQNe provide ample ioniz- 
ing photons per baryon to re-ionize the universe and provide the 
optical depth in the CMB of ~ 0.1 detected by WMAP. Thus ac- 
cording to our model, the end of the epoch of re-ionization should 
also coincide with the epoch of formation of Pop II low-mass stars. 
A notion which could be verified observationally. 
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Figure 12. Relative abundances [X/Fe] of sub-Fe spallation products (elements identified by their names) versus mass number A from simulations with ^dclay = 
3.5 days which resulted in [Fe/H] = -4.8. The CN enhanced nitrogen is shown by the filled diamond. Left and right panel are for QN neutrons with Eqn = 5 
GeV and 10 GeV, respectively. SDSS J102915-l-172927's measured abundances are also shown, as plus (-I-) signs with the corresponding uncertainties, for a 
comparison (plus signs with no error bars are upper limits). We adopted measured abundances given in Caffau et al. (201 1 ; see their Table 1) with an assumed 
upper limit for oxygen of [O/Fe] ~ 0.6. Best fits were found for initial abundances in mass of Mni.sn = 0. IMq, Ms,sn = ^Si.SN = /^Mg.sN = /Wnc.sn = O.OIMq, 
A^o.SN = O.O4M0 and Mc.sN = 0.0.02Mo. The mass of the cloud swept by a dsQN was kept at Ms„ = IO^'Mq. 



6 CONCLUSION 

In this paper we made a case for QNe (specifically dsQNe) as 
sources of ^Li in the pre-Galactic era. Lithium is produced by spal- 
lation on the PopIII-SNII ejecta material. For dsQNe with ?delay > 
Pouter ~ 8.3 days, an A(Li) ~ 2.2 ^Li plateau results from spallation 
in the Oxygen layer matching observations of Spitec&Spite (1982). 
We argue for a synergy between our model and the astration model 
of Piau et al. (2006). In the astration-dsQN scenario, the BBN is 
first processed/destroyed in the 20-40Mo Pop III QN progenitor 
before it is restored to Spite values following spallation in oxygen 
layer of the PopIII-SNII ejecta. In our model, instabilities in the 
pristine cloud swept-up by the SN+QN (i.e. dsQN) ejecta should 
lead to the formation of low-mass stars with abundances reminis- 
cent of those measured in old metal-poor dwarfs and in turn-off 
stars of globular clusters. According to our model, the most iron- 
poor stars likely formed from gas clouds swept up and enriched by 
dsQNe with very short delays. This effectively dissociates age from 
metallicity and thus calls for a paradigm shift in the current wisdom 
which seems to associate the iron content in stars increases with the 
cosmic age. 

In today's universe, QNe should be associated with mas- 
sive star (2O-4OM0) regions. At low redshift, dsQNe would 
manifest themselves as superluminous SNe (Guyed et al. 2002; 
Leahy&Ouyed 2008; Ouyed et al. 2010) with a distinct "double- 
hump" light-curve - the first hump correspond to the core-collapse 
SN explosion proper and the second hump to the re-energization 
of the SN shell by the QN (see Ouyed et al. 2009b; Ouyed&Leahy 
2012). Other imprints of dsQNe in today's universe have been sug- 
gested in the context of Ni-poor, Ti-rich SNe such as Gas A and in 
the context of SNe with peculiar chemical composition in general 
Guyed et al. (2011). Interestingly, we find a similar range in fjelay 
(2 days < fdelay < 20 days) when fitting SLSNe light-curves (e.g. 
Guyed et al. 2012; Kostka et al. 2012) and the chemical composi- 
tion of some very peculiar SNe (Ouyed et al. 201 1). This might be 
an indication (or plausibly a self-consistency check) on the univer- 
sal nature of the physics underlying the QN explosion and of the 
mass of their progenitors. This range in /delay combined with the 
estimate of the transient gravitational wave signal from QNe (Staff 



et al. 2012) should lead to distinct dsQNe signatures in future Grav- 
itational Waves detectors. 

The existence of deconfined quark matter in the superdense 
interior of neutron stars is a key question that has drawn consider- 
able attention over the past few decades. Here we have shown that 
the existence of the Spite plateau is an indication that quark matter 
does appear in the core of massive neutron stars and that it triggered 
the biggest explosions in the universe, namely the QNe. These QNe 
may also have played a role in the early universe by helping in its 
re-ionization (Ouyed et al. 2009a) and by triggering the formation 
of very low-metal Pop II stars (Ouyed 2013 and this work). Other 
implications to Cosmology have been discussed in Ouyed&Staff 
(2013). 
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